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BOUNDARY CONDITIONS FOR THE SINGLE-FACTOR TERM 
STRUCTURE EQUATION 

By Erik Ekstrom^'^ and Johan Tysk^ 

Uppsala University 

We study the term structure equation for single-factor models 
that predict nonnegative short rates. In particular, we show that the 
price of a bond or a bond option is the unique classical solution to a 
parabolic differential equation with a certain boundary behavior for 
vanishing values of the short rate. If the boundary is attainable then 
this boundary behavior serves as a boundary condition and guaran- 
tees uniqueness of solutions. On the other hand, if the boundary is 
nonattainable then the boundary behavior is not needed to guaran- 
tee uniqueness but it is nevertheless very useful, for instance, from a 
numerical perspective. 

1. Introduction. When calculating prices of different interest rate deriva- 
tives, such as bonds and bond options, stochastic methods seem to be more 
commonly used than PDE methods; compare for instance [5]. This is in 
contrast to the case of stock option pricing where PDE methods are used 
extensively, in particular for low dimensional problems. We believe that one 
possible explanation for this phenomenon is that the correspondence be- 
tween the risk neutral valuation approach and the pricing equation (hence- 
forth referred to as the term structure equation) with appropriate boundary 
conditions is not fully developed. 

For the Black-Scholes equation, the boundary condition is of Dirichlet 
type which corresponds to the underlying asset being absorbed if reaching 
zero; compare [13]. In contrast, for most interest rate models this is not 
the case since the short rate typically would not stay zero if the value zero 
is reached. Consequently, it is not clear what boundary conditions should 
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be specified for the term structure equation. In fact, the recent monograph 
[6] draws attention to this issue in a section entitled "The thorny issue of 
boundary conditions." Moreover, in [12], two different solutions to the term 
structure equation are presented in the case of the CIR-model. They are 
both bounded and with the same terminal condition, but naturally they ex- 
hibit different boundary behavior for vanishing interest rates. The authors of 
that paper take the view that these solutions represent alternative possible 
prices. We on the other hand regard only the solution given by the stochas- 
tic representation as the price, and the purpose of the present paper is to 
identify the boundary condition that this stochastic representation satisfies. 

We consider the classical case of a single- factor model predicting non- 
negative values of the short rate. More precisely, the rate X{t) is modeled 
directly under the pricing measure as 

dX{t) = l3{X{t),t) dt + a{X{t),t) dW, 

where is a Brownian motion and o"(0,t) = and /3(0,t) > 0. As indicated 
above, the option price u corresponding to a payoff function g is given using 
risk neutral valuation by 

u{x,t) = E^4e~f"''^'^^'g{X{T))]. 

Note that if the payoff g =1, then bond prices are obtained. Also note that 
the set-up covers the case of bond options. The corresponding term structure 
equation is 

Ut{x, t) + |cT^(x, t)Uxxix, t) + (3{x, t)Ux{x, t) = Xu{x , t) 

with terminal condition u{x,T) = g{x). If the price u is twice continuously 
differentiable up to and including the boundary x = 0, then plugging in x = 
in the equation would give the boundary behavior 

(1) nt(0,t) + /3(0,t)n^.(0,t) = 0. 

Even though there is extensive literature on equations with degenerating 
coefficients, compare the classical reference [14], the C^-regularity of u at 
the boundary is not available in the generality that is needed here. In fact, 
one of the solutions in the example in [12] referred to above is bounded and 
continuous, but fails to be up to the boundary. (This solution, however, 
is not the one given by stochastic representation.) 

In the present paper, sufficient regularity of the option price u for (1) to 
hold is established using the Girsanov theorem and scaling arguments; see 
Sections 3 and 4. Let us emphasize that (1) is the correct boundary behavior 
of the option price regardless if the boundary is hit with positive probability 
or not. If the boundary can be reached with positive probability, then this 
boundary behavior serves as a boundary condition for the term structure 
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equation and guarantees uniqueness. On the other hand, if the boundary 
is reached with probabihty zero, equation (1) is not needed to identify the 
solution given by the stochastic representation, and the term "boundary 
condition" is perhaps misleading. However, it is still valid and certainly 
useful, for instance, from a numerical perspective. Indeed, the results of 
the present paper have already been implemented numerically in [7]. For 
simplicity of the exposition, we will refer to (1) as the boundary condition 
regardless if the boundary can be reached or not. 

In Section 2 the assumptions on the model and our main result. Theo- 
rem 2.3, are presented. In Sections 3, 4 and 5, we establish regularity prop- 
erties of the value function, and use them to prove Theorem 2.3. Finally, in 
Section 6 we also provide the link between the stochastic problem and the 
term structure equation for models defined on the whole real line. 

2. Assumptions and the main result. When studying the term structure 
equation on the positive real axis, we consider models which are specified 
so that the short rate automatically stays nonnegative, that is, there is no 
need to impose any boundary behavior of the underlying diffusion process. 
Throughout Sections 2-5 we work under the following hypothesis: 

Hypothesis 2.1. The drift f3 G C([0,oo) x [0,r]) is continuously differ- 
entiable in x with bounded derivative, and /3(0,t) > for all t. The volatility 
a G C([0,oo) X [0,T]) is such that a{x,t) := ^a'^{x,t) is continuously differ- 
entiable in x with a Holder continuous derivative, and a{x,t) = if and only 
if x = 0. The functions (3, a and Ox are all of, at most, linear growth: 

(2) \p{x,t)\ + \a{x,t)\ + \ax{x,t)\<C{l + x) 

for all X and t. The payoff function g: [0,oo) — )■ [0, oo) is continuously dif- 
ferentiable with both g and g' bounded. 

Let W he a standard Brownian motion on a filtered probability space 
{Q,T,(J^t)t>o,P)- Since a is continuously differentiable, a is locally Holder 
(1/2) in X. It follows that there exists a unique strong solution X{t) to 

(3) dX{t) = P{X{t),t) dt + (T{X{t),t) dW 

for any initial point x > 0; compare Section IX. 3 in [15]. Moreover, it follows 
from monotonicity results for stochastic differential equations with respect 
to the drift coefficient (e.g.. Theorem IX. 3. 7 in [15]), that X remains nonneg- 
ative at all times. Indeed, if (3 is replaced with /? AO, then the corresponding 
solution to (3) is absorbed at zero, so X is nonnegative if /3(0,t) > 0. The 
option price u: [0,oo) x [0, T] — )• [0,oo) corresponding to a payoff function 
g : [0, oo) — )• [0, oo) is given by 

(4) uix,t) = Ex,t[e-^"''^'^'''giXiT))], 
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where the indices indicate that X{t) = x. As described in the Introduction, 
the corresponding term structure equation is given by 

(5) ut{x,t) + ^cr^(a;,t)n^.^.(x,t) + (3{x,t)ux{x,t) =xu{x,t) 
for {x,t) £ {0,00) X [0,r), with terminal condition 

(6) u{x,T) = g{x). 

Moreover, by formally inserting x = in the equation we get the boundary 
condition 

(7) nt(0,t)+/3(0,tK(0,t) = 

for all t G [0,T), since o"(0,t) = by assumption. One of the main efforts in 
this paper is to show that the option price u is continuously differentiable up 
to the boundary x = 0, and that it indeed satisfies the boundary condition 

(7) in the classical sense. 

Definition 2.2. A classical solution to the term structure equation is 
a function v£C{[0,oo) x [0, T]) n C^([0, 00) x [0, T)) n C2'i((0, 00) x [0,T)) 
which satisfies (5), (6) and (7). 

Our main result in this article is the following: 

Theorem 2.3. In addition to Hypothesis 2.1, also assume that Assump- 
tion 3.1 below holds. The option price u as given by (4) is then the unique 
bounded classical solution to the term structure equation. 

Example. Classical short rate models such as the Cox-Ingersoll-Ross 
model 

(8) dX{t) = {a- bX{t)) dt + a^JX{t) dW, 
and the Dothan model 

(9) dX{t) = aX{t)dt + aX{t)dW, 

have boundary conditions at x = that are immediate to write down. These 
conditions are 

ut + aux = and ut = 0, 

respectively. We note that the boundary condition = for the Dothan 
model means that u is constant along the boundary, that is, ?i(0,t) = g{Q). 
This is the same type of boundary condition that appears for options on 
stocks in [13], which can be explained by the fact that the Dothan model 
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is a geometric Brownian motion. Theorem 2.3 also covers the Hull-White 
model 

(10) dX{t) = {a{t) - b{t)X{t)) dt + a{t)y^X{t)dW 

(which is a time-dependent generalization of the Cox-Ingersoll-Ross model), 
and models of, for example, the form 

(11) dX{t) = {b-aX{t))dt + aX^{t)dW, 7 G (1/2,1], 
which also would be natural to consider for bond pricing. 

Remark. It seems that many of the classical models for the short rate 
are proposed for their analytical tractability. In particular, if the drift /3 
and the diffusion coefficient cr^ are affine, then the model admits an affine 
term structure. It is easy to check that known explicit formulas for bond 
prices and bond options satisfy the boundary condition (7). In particular, 
for models admitting an affine term structure, it is a consequence of the 
associated Riccati equations (see [4], equation 22.25) that these boundary 
conditions are fulfilled. 

Remark. The assumption that g is continuously differentiable is sat- 
isfied for bonds, but not in general for bond options. However, using the 
Markov property. Theorem 2.3 readily extends to bounded Lipschitz payoffs 
provided one can show that the corresponding option price x u{x,T — e) 
is continuously differentiable on (0, 00) for any e > 0. The regularizing effect 
of parabolic equations guarantees continuous differentiability on (0,oo), so 
the main difficulty is to show that x Ux{x,T — e) is continuous also at 
0. If the model is convexity preserving, this is easily done in certain cases 
including, for example, call options written on bond prices. (Note that the 
corresponding payoff function g is bounded since bond prices are bounded.) 
For details on which short rate models are convexity preserving, see [8]. To 
our knowledge, all models used in practice belong to this class. 

One should note that the differentiability of the option price up to the 
boundary x = is not valid without some Lipschitz bound of g at 0. To see 
this, consider the contract function g{x) = e~^^. Then, with (3{x,t) = ^x 
and a{x,t) = y^x, it is straightforward using the Ito formula to show that 
the process 

y(,) = e-//^M'^'-5(X(s)) 

is a martingale. Consequently, the option price u is given by u{x,t) = g{x) 
for all t, which fails to be a classical solution to the term structure equation 
since it is not differentiable at x = Q. One might argue, though, that the 
boundary condition ut{0,t) + f3{0,t)ux{0,t) = is satisfied in a weak sense. 
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The proof of Theorem 2.3 is carried out in several steps. 

Proof of uniqueness. Let and v"^ be two bounded classical solu- 
tions to the term structure equation, and define 

v{x,t) = v^{x,T - t) - v^{x,T - t). 

Then v{x,t) is a bounded solution to 

{vt = ^a^Vxx + - XV, 
v{x,0) = 0, 
vt{0,t)=(3{0,t)vx{0,t). 

Now consider the function 

h{x,t) = {l + x)e^\ 

where M is a positive constant. For M large enough, depending on /3(0,t) 
and the growth rate of /3, /i is a super-solution to (12) which tends to infinity 
at spatial infinity. Thus, according to the maximum principle, the function 
V is bounded above by eh and below by —eh for any e > 0. It follows that 
v = 0, which demonstrates uniqueness of bounded classical solutions to the 
term structure equation. □ 

Proof of continuity. To show that u is continuous, denote by X^'* 
the solution to (3) with initial condition X^'*(t) =x. Let {x,t) and (y,r) be 
two points in [0,oo) x [0,T]. Then, if r < t, we have 

\uiy,r)-u{x,t)\<E[e~f"'''-''^'^'''\g{Xy'^{T))-g{X-'\Tm 
(13) <E[\g{Xy^'{T))-g{X-'\Tm 



+ C E[\Xy^''{s)-X''^\s)\]ds 
+ C f E[Xy'' {s)]ds 

J r 



for some constant C, where we have used that g is bounded. A similar 
expression can be derived if r > t. It follows from Remark 1 in Section 8, 
Chapter 2 in [11] that X^''''(t) — )■ x in as (y,r) — )■ (x,t). Therefore, from 
Theorem 2.1 in [2] we have 



E 



sup {Xy''{s)-X'''\s)f 

t<s<T 







as (y, r) — )• (x, t). (Theorem 2.1 in [2] also holds in the case of random starting 
points.) Since g is assumed continuous and bounded, all three terms on the 
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right-hand side of (13) tend to as (y,r) — )• (x,t). Thus u is continuous on 
[0,oo) X [0,T]. □ 

Proof that u£ C^'-'^((0,oo) x [0,T)) and satisfies (5). For a given 
point (x,t) G (0,oo) x [0,r), let 

R={xi,X2) X [tiM) ^ (0,oo) X [0,r) 

be a rectangle which contains {x,t), where xi > 0. Since u is continuous, it 
follows from standard parabolic theory, see [9], that there exists a unique 
solution U G C^'^(i?) to the boundary value problem 

iUt + ^a^U^^ + ^U^-xU = 0, mR, 
\U = u, on dpR, 

where dpR = {[xi,X2] x {t2}) U ({xi, 3:2} x [ti,t2]) is the parabolic boundary 
of R. From Ito's formula, the process 

Z{s) = e- '^'U{X'''\s),s) 
is a martingale on the time interval [t,TR], where 

TR = mi{s>t:X'''\s)(^R} 
is the first exit time from the rectangle R. Therefore, 

U{x, t) = E[e- f^^"" ''^u{X-'\tr),tr)] = u{x, t), 

where the second equality follows from the strong Markov property. Con- 
sequently, u G C^'^((0,oo) X [0,T)). Since u = U on R, we also see that u 
satisfies (5). □ 

It remains to show that u is continuously differentiable up to the spatial 
boundary x = 0, and that it satisfies the boundary condition (7). This is 
done in Sections 3-5. 

3. Continuity of the first spatial derivative. In this section we investi- 
gate regularity of the spatial derivative Ux at the boundary x = 0. To do 
this we study the stochastic representation of the terminal value problem 
obtained by formally differentiating the term structure equation. We show 
that this stochastic representation indeed is the derivative of u and that it 
is continuous. 

Recall that Ox is assumed to be continuous on [0,cx3) x [0,T], where 
a{x,t) = ^o"^(x,t). Let the process Y be modeled by the stochastic differen- 
tial equation 

(14) dY{t) = [ox + ^){Y{t),t) dt + (T{Y{t),t) dW. 

Rather than specifying precise conditions under which (14) has a unique 
solution, we simply assume what we need. 
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Assumption 3.1. The coefficients a and (3 are such that, path-wise, 
uniqueness holds for equation (14)- 

Remark. Note that Assumption 3.1 holds for example if a is twice 
continuously differentiable in space, since then the drift + /3 is locally 
Lipschitz continuous. Moreover, if a and f3 are time-independent, then it 
follows from [1, 3] and Section IX. 3 in [15] that Assumption 3.1 automatically 
holds. Thus the Cox-Ingersoll-Ross model (8), the Dothan model (9), the 
Hull-White model (10) and the model (11) all satisfy Assumption 3.1. 

Also note that since a{0,t) = 0, we have ax{0,t) > 0. Thus Y remains 
nonnegative since it has the same volatility as X but a larger drift at 0. 

Next, define the function v by 

vix, t)=E g'(y(T)) expj j%x{Y{s),s) - Y{s) ds 

(15) 

-E 



^ exp|^ f3,{Y{r),r)-Y{r)dr^u{Y{s),s)ds 



where Y is the solution to (14) with initial condition Y{t) = x. 

If the term structure equation (5) is formally differentiated with respect 
to X, then the derivative Ux satisfies 

{Ux)t + Ol{Ux)xx + {Oix + l3){Ux)x + ii3x - x)Ux - U = 

with terminal condition Ux{x,T) = g'{x). The function v defined in (15) is 
the corresponding stochastic representation. In Theorem 3.4 below we show 
that V indeed equals the spatial derivative of u. 

Proposition 3.2. The function v{x,t) is continuous on [0,oo) x [0,T]. 

Proof. The result follows along the same lines as the continuity of u 
above. Indeed, let {xn,tn) converge to {x,t), where tn < t, and let Y and y" 
be defined by 

f dY{s) = {ox + /3)(y(s), s) ds + a{Y{s),s) dW, 
\Y{t)=x 



and 



respectively. Also define 



(iy"(s) = (ax + /3)(y"(s), s) ds + cj(y"(s), s) dW, 

y (^n) — ^r, 



I{s) :=exp|^'/3. 



{Y{u),u)-Y{u) du 
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and 



/"(s):=exp<{ / /3^(y"(u),n) -y"('u)dn 



Then 



\v{xn,tn)-vix,t)\ < E[\P {T)g' [Y^iT)) - I{T)g' {Y {T))\] 

+ /' E[\P{s)u{Y^{s),s)-I{s)u{Y{s),s)\]ds 



+ E[r{s)u{Y'\s),s)]ds. 

J tn. 



The first term and the integrand in the second term are similar to the type 
of terms treated when proving the continuity of u. Moreover, the integrand 
of the third term is bounded. Thus it follows from bounded convergence that 
V is continuous. □ 



We also need a continuity result in the volatility parameter. To formulate 
it, let {(t"(j;, t)}^^ be a sequence of functions satisfying Hypothesis 2.1 
uniformly in n, that is, with the same constant C in the bound (2). Moreover, 
assume that a'"'{x,t) converges to a{x,t) and a" converges to Ox uniformly 
on compacts as n — t- oo, where a" = ^(cj")^. Let n" and v"" be defined as u 
and V but using the volatility function o"" instead of a. More explicitly. 



and 

v''{x,t)=E 
(16) 



g'iY^{T))exp{ / l3xiY^{s),s) -Y^s) ds 



E 



^ exp I (y" (r ) , r ) - y" (r ) [> n" (y" (s) , s) 



where and y" satisfy 



and 



dX"(s) = /3(X"(s), s) ds + cj"(X"(s), s) dW{s) 



(iy"(s) = « + /3)(y"(s), s) ds + a"(y"(s), s) dWis), 
y"(t) =rE, 



respectively. 



10 E. EKSTROM AND J. TYSK 

Proposition 3.3. The functions u andv are continuous in the volatility 
parameter. More precisely, u^{x, t) — )■ u{x, t) and v^{x, t) — )• v{x, t) as n —)• oo 
for any fixed point {x,t) £ [0,oo) x [0,T]. 

Proof. It follows from Theorem 2.5 in [2] that 

= 0. 



lim E 

n— >oo 



sup (X(s)-X"(s))^ 

-selt,T] 



Therefore, 

\u''{x,t)-u{x,t)\ < ^[|e-/t^^"('')'^^-e-ii^^(^)'^^|5(X'^(T))] 



+ ^[e-/"^(^)^^|9(X"(r))-g(X(r))|] 

< C j%[\X{s) - X^{s)\] ds + EMX-{T)) - g{X{T))\] 
^0 

as n — 7- oo. Thus u is continuous in the volatility function. 

The continuity of v in the volatility function is similar. Indeed, let 

I{s):=exp!^j^ /3,(y(r),r)-y(r)dr| 

and 

/" (s) := exp I (y" (r ) , r) - (r) dr | . 

Then 

\v-ix,t)-vix,t)\<E[\P{T)g\Y-iT))-I{T)g'{Y{T))\] 

+ r E[\ris)u^{Y^is),s)-Iis)u^iYis),s)\]ds 



+ E[Iis)\u^iYis),s)-uiYis),s)\]ds. 

The first term and the integrand of the third term are similar to the terms 
appearing in the first part of the proof. The integrand of the second term 
can be dealt with using the fact that each it" is Lipschitz continuous in x, 
uniformly in n since the Lipschitz property is inherited by the value function. 
Thus all terms tend to zero as n — >• oo, so f is continuous in the volatility. 
□ 



Theorem 3.4. We have Ux{x,t) = v{x,t) on [0,oo) x [0,r]. Conse- 
quently, Ux is continuous on [0,oo) x [0,T]. 
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Proof. It suffices to prove Ux{x,0) = v{x,0). We first assume that a is 
continuously differentiable in x with a bounded derivative. It then follows 
from Section 5.5 in [10] or Section 8 in [11] that the derivative 



dXjt) 
dx 



of X{t) = X^'^{t) with respect to the initial point x exists and is continuous, 
and it satisfies 

/ dm = m^x{x{t),t) dt + t) dw{t), 

U(o) = i. 



Moreover, 



Ux{x,0) = E 



g'{X{TmT)exp 



X{s)ds 



(17) 



-E 
h-h. 



<7(X(r))exp|-^ X{s)dsjl^ C{s)ds 



We claim that /j = Jj, i = 1,2, where 



Ji = E 



g'(y(r))exp / /3x{Y{s),s)-Y{s)ds 



and 



J2=E 



exp 



/3^(y(r), r) - y(r) dr \u{Y{s),s) ds 



compare (15) above. Here Y is defined as in (14) with initial condition 

y(o) = X. 

To show that /i = Ji, define a new measure Q on J^x by dQ = M(T) dP, 
where the process M is defined by 



(18) M(t)=e(t)exp|- /3x{Y{s))ds 

By Ito's formula, 

dM{t) = M{t)cTx{X{t)) dW{t), 

so M is a martingale since Gx is bounded. In particular, E[M{T)] = 1, so Q 
is a probability measure. From Girsanov's theorem it follows that 



W{t) = W{t)- f ax{X{s))ds 
Jo 
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is a Q-Brownian motion, and 

dX = {aa^ + P)iX{t),t) dt + a{X{t),t) dW. 

Here aa^ = Ox, so by weak uniqueness, the Q-law of X is the same as the 
law of Y under P. Consequently, 



E 



g'(x(r))r(r)exp 



X(s)ds 



eQ 



g'iX(T))exp\ / f3xiX{s), s) - X{s) ds 



Ji. 



To prove l2 = J2, note that 



■ E 



5(X(r))exp|-^ X{s)dsjl^ as)ds 



E 



exp 



X E 



X{r)dr }^{s) 



<7(X(r))exp 



X[r) dr 



X{r)dr\i{s)u{Xs,s) 



ds 



ds 



by the Markov property. Define a new measure Q = Qs on J^s by 

dQ = M{s) dP, 
where M is defined as in (18). Girsanov's theorem yields 



E 



exp 



E^ 



exp 



X{r)dr'j^{s)u{Xs,s) 

(3x{X{r),r)-Xir)dr}u{Xs,s) 



E 



exp 



(3x{Y{r),r)-Y{r)dr}u{Ys,s 



Consequently, I2 = J2, which finishes the proof in the case of continuously 
differentiable a. 

The general case follows by approximation. Let o"", and be as 
described before Proposition 3.3, with each o"" being continuously differ- 
entiable in X with bounded derivative. From above, we then know that 
v^{x,t) = vJ^{x,t) at all points. Moreover, by Proposition 3.3, v"'{x,t) — )• 
v{x,t) point-wise as n — )• 00. 
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On the other hand, since converges to u point-wise and is uniformly 
bounded, it foUows from standard parabohc theory that also converges 
to Ux point-wise for all points {x,t) with x > 0. Consequently, v = Ux on 
(0,00) X [0,T]. Since v is continuous on [0,oo) x [0,T] by Proposition 3.2, it 
is easy to check that ^^(O,^) exists and that we have v = Ux everywhere on 
[0,00) X [0,T]. The continuity of Ux thus follows. □ 

4. An estimate of the second spatial derivative. Since the function v 
defined in (15) is continuous, it follows that [by a similar argument as in the 
proof that u satisfies (5)] it indeed solves the differentiated equation 



on (0,00) X [0,T). In this section we use interior estimates to show that 
avx — 7- as a; — ?• 0. Since v = Ux by Theorem 3.4, this shows that the term 
auxx in (5) approaches zero close to the boundary. 

Proposition 4.1. The function v = Ux satisfies 



for any to- Consequently, lim(^x,t)^{o,to) <^{x,t)uxx{x,t) =0. 

Proof. Let {{xn,tn)}^^i ^ (0,oo) x [0,T) be a sequence of points con- 
verging to (0,to)5 where to S [0,T). Define new coordinates by letting 
y = kx and s = k{t — to), where k is specified more precisely below. Then 
the function w defined by 



Vt = aVxx + {otx + P)Vx + (A 



'X 



x)v — u 



lim a{x,t)vx{x,t) = 

(x,f)-s>(0,to) 



w{y,s) =v{x,t) 



satisfies 



(19) 



■Ws = OiWyy + f3Wy +JW + h, 



where 



a 






and 
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Now consider a region 71 = 7^^ which contains the point {xn,tn), and such 
that 

(20) l<a{x,t)k<2 

in TZ. Since ax{x,t) is continuous up to the boundary, the region TZ in 
(y, s)-coordinates does not collapse as n— )-oo, but it can rather be chosen 
to consist of a rectangle of fixed size; the location of the rectangle is not 
necessarily fixed though. In this rectangle, the coefficients of the equation 
(19) satisfy 

l<a{y,s)<2, 
my,s)\<c, 
hiy,s)\<C 

and 

\Ky,s)\<C/k 

for some constant C which is independent of n. Since w{y,s) =v{x,t) we 
have that w converges to the constant u(0,io) = Ux{0,to) uniformly on TZ as 
n — )• oo. By interior Schauder estimates, Wy tends to as n — )• oo. Since 

a{x,t)vxix,t) = a{y,s)wy{y,s), 

and since a{y, s) is bounded on TZ, the conclusion follows. □ 



5. The time derivative at the boundary. It follows from Proposition 4.1 
and (5) that 

(21) lim ut{x,t) + (3{0,to)uxiO,to) = 

for any to £ [0)^)- In this section we show that the boundary condition (7) 
also holds at the boundary, that is, not merely in the limit. 

Proposition 5.1. The function ut{x,t) + P{x,t)ux{x,t) defines a con- 
tinuous function on [0,oo) x [0,T). Moreover, it vanishes for x = 0. 

Remark. Note that Proposition 5.1 finishes the proof of Theorem 2.3. 

Proof of Proposition 5.1. In view of (21) above, it suffices to show 
that ut exists at the boundary and that it equals —(3ux- To do this, fix a 
point on the boundary with coordinates (0,to)- For notational simplicity we 
assume that to = 0. The time (left) derivative ut at the boundary is defined 
by 

(22) u^(o,0)= lim /c( ti(0,0) -uf 0,-i 

k~^oo \ \ k 
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provided the limit exists. To determine ut{0,0), we let be defined by 

dX'' =(3{X''{t),t)dt + a{X''{t),t)dW, 
X''{-l/k) = 0. 

However, instead of considering the process X with different starting times, 
we perform a change of variables so that the starting time is independent of 
k. We thus introduce the process Y^{s) by 



Y\s) = kX^ J-]. 



With respect to the time variable s, the dynamics of has the form 
(23) 



[yfc(-l) = 0. 



where W^{s) denotes some Brownian motion. By the Markov property, 
u{0, -l/k) = E[e~^-^/^^^'^'^'^%{X''{0),0)] 



■ E 



Hence, 



nt(0,0) = lim fc-Bo,-! 



: lim Eq^^i 



k( u{0,0)-u( -Y''{0),0 



where the second equality follows using the inequality e ^ — 1 > —x since 



ku{ ^y'=(0),0 )|e- 



k 



Y^{s)ds 



as fc — )• oo. Now, define the process Y by 



fdY = /3(0, 0) ds + y/2a^ {0,0)Y dW, 

\y{-i) = o, 

and redefine Y^ as in (23) above but using the same Brownian motion W 
(this does not change the law of Y''). Since 
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and 



y s 
k' k 



\/2a.(0,0)y 



uniformly on compacts as /c — )• oo (here we used the assumption that a is 
continuously differentiable in space), it follows from [2] that Y^{0) — )• 1^(0) in 
as A; — 7- oo. From Theorem 3.4 above we know that u is differentiable in x, 
so A;(n(0,0) — ii(|,0)) converges to —Ux{0,0)y. By dominated convergence, 
we have 



kE 



u{0,0)-u{ ^Y{0),0 



-u,{0,0)E[Y{0)] = -/3(0,0)n,(0,0) 



as fc — 7- oo. Moreover, the Lipschitz property of u yields that 



kE 



n(^y(0),0 



n(iyno),o 



<c^[|y(o) -y'=(o)|] ^0 



as k 



oo. It follows that 

ut(0,0)+n^(0,0)/3(0,0)=0. 

As to = was chosen only for notational convenience, we have that 

'Ut(0,t) + /3(0,tK(0,t) = 

for any t. To be precise, we have shown the result above only for the left 
t-derivative. However, this left t-derivative is continuous by the equation 
above, so it follows from a simple calculus lemma that in fact u is differen- 
tiable in time, thus finishing our proof. □ 



6. Models allowing negative interest rates. For models in which the 
short rate can fall below zero with positive probability, the connection be- 
tween the option price, given by a risk-neutral expected value, and the term 
structure equation is more straightforward than for models with nonnega- 
tive rates. Nevertheless, we have not been able to find a precise reference for 
this case, so for completeness we provide such a result in this section. 

The assumptions needed on the drift and volatility are presented below, 
where = max(x,0). These assumptions now replace those of Hypothe- 
sis 2.1. To our knowledge, all models used in practice that allow negative 
interest rates satisfy these requirements. For example, the Vasicek model, in 
which 

dX{t) = {a- bX{t)) dt + a dW, 
where a, b and a are positive constants, is covered. 
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Hypothesis 6.1. The drift (3 G C(M x [0,r]) and the volatility a € 
X [0,T]) are both Lipschitz continuous in x. Moreover, 



0< |cr(2;,t)| <C(l + a;+) 

and 

(24) \l3{x,t)\<C{l + \x\) 

for some positive constant C . 

Remark. The assumption that a is strictly positive is not a strong 
assumption. Indeed, let us for simplicity consider a time-homogeneous model 
dX{t) = /3{X{t))dt + a{X{t)) dwlt) with a{a) = for some a G M. If X(0) > 
a and /3(a) > 0, then the process X cannot take values smaller than a, so we 
are essentially in the situation handled in Sections 2-5 (but with the point 
replaced with a). If 6(a) < 0, then X can take values below a, but if this 
happens then the process will stay below a forever, and we are then again 
back in the previous situation. 

The bound on the volatility for negative rates guarantees that bond prices 
are finite. For models in which a grows faster than y^|x| for negative rates, 
bond prices can be infinite; compare Theorem 4.1 in [16]. 

For a given continuous payoff function g:R— t- [0,oo), define the corre- 
sponding option price n : M x [0,T] by 

/•T 

u{x,t) = E^^t 



expj-^ X{s)ds^g{X{T)) 



where 



dX{s) = P{X{s),s) ds + a{X{s),s) dW{s), 
X{t)=x. 

We require that the payoff function is bounded for positive interest rates 
and of, at most, exponential growth for negative rates, that is, 

(25) 0< c/(x) <Kmax{l,e"-^'^} 

for some positive constant K. The corresponding term structure equation is 
given by 

ut{x,t) + ^a'^{x,t)uxx{x,t) + f3{x,t)ux{x,t) = xu{x,t) 

on M X [0,T), with terminal condition u{x,T) = g{x). By a classical solution 
to this equation we mean a solution which is continuous up to the boundary 
t = T and with all derivatives appearing in the equation being continuous 
functions on the set t <T. The following is our main result about the term 
structure equation on the whole real line. 



18 



E. EKSTROM AND J. TYSK 



Theorem 6.2. Assume Hypothesis 6.1 and the hound (25). Then the 
option price u{x,t) satisfies 

(26) u{x,t)<K'max{l,e-^'''} 

for some constant K' . Moreover, u{x,t) is the unique classical solution to the 
term structure equation satisfying this growth assumption for some constant 
K'. 

Remark. Note that the bound (25) is natural for models on the whole 
real line. In fact, even if g was bounded, the option price u would be of expo- 
nential growth for negative rates. Also note that, for example, call options 
on a bond are covered by Theorem 6.2. Indeed, the payoff of a bond call 
option with maturity Ti is given by g{x) = {u{x,Ti) — K)^ , where u is the 
price of a bond maturing at T2 > Ti. Since u satisfies (26) by Theorem 6.2, 
the payoff g satisfies (25). Bond put options are trivially covered since they 
have bounded payoff functions. 

Proof of Theorem 6.2. The bound (26) follows from Corollary 3.3 
in [8]. To prove uniqueness of solutions, assume that w is a solution to 
the term structure equation with boundary value g = such that \v\ < 
i^'max{l,e~^ ^} for some constant K'. Let 

h{x,t) = e^'^^-'\e-f^'>+x) 

for some large constant M. Here 

pC(T-t) 1 

f{t)= ' ^-^ Ke^^^-^ 
where C is the constant appearing in (24) and K > K' . Then the set 

{(x,t) gM X [{),T]:eh{x,t) <v{x,t)} 

is bounded, and 

ht + \o''^hxx + — x/i < 

at all points provided M is chosen large enough. Standard methods used 
to prove the maximum principle yield that v = Q aX all points. Thus we 
have uniqueness of solutions to the term structure equation in the class of 
functions satisfying (26). 

To show that u is a classical solution to the term structure equation, 
we carry out an approximation argument. We consider the term structure 
equation but with the discount factor x replaced by bounded functions that 
agree with x inside some large compact sets. We also replace the payoff 
function g with functions of, at most, polynomial growth that agree with g 
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on large compact sets. The corresponding equation is in the standard class 
and its stochastic solution is known to be continuous and hence is a classical 
solution. Now let the functions approximating the discount factor grow up 
to X for large positive x and then decrease down to x for large negative x. 
By the monotone convergence theorem, the corresponding stochastic solu- 
tions converge to the stochastic solution above denoted u. Interior Schauder 
estimates yield interior regularity of the limiting solution and continuity at 
the boundary is established using the maximum principle. □ 
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